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Abstract: 

We investigate an effective model for the finite temperature symmetry 
restoration phase transition of the electroweak theory It is obtained by di- 
mensional reduction of the 3 + 1 dimensional full theory and by subsequent 
integration over all static gauge degrees of freedom. The resulting theory cor- 
responds to a 3-dimensional 0(4) ferromagnet containing cubic and quartic 
terms of the field in its potential function. Possible nonperturbative effects 
of a magnetic screening mass are parametrically included in the potential. 
We analyse the theory using mean field and numerical Monte Carlo (MC) 
simulation methods. At the value of the physical Higgs mass, m H = 37 GeV, 
considered in the present investigation, we find a discontinuous symmetry 
restoring phase transition. We determine the critical temperature, order pa- 
rameter jump, interface tension and latent heat characteristics of the transi- 
tion. The Monte Carlo results indicate a somewhat weaker first order phase 
transition as compared to the mean field treatment, demonstrating that non- 
perturbative fluctuations of the Higgs field are relevant. This effect is espe- 
cially important for the interface tension. Any observation of hard first order 
transition could result only from non-perturbative effects related to the gauge 
degrees of freedom. 
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1 Introduction 



The tiny matter-antimatter asymmetry observed on cosmic scales requires 
explanation in the framework of standard cosmology. Baryon number vi- 
olating processes are known to occur in the electro-weak theory and are, 
in fact, expected to occur frequently at high temperatures in equilibrium. 
The occurence of a non-vanishing matter-antimatter asymmetry, however, 
is expected to be possible only, if the universe evolved through some non- 
equilibrium stages. If the electroweak phase transition was of first order 
(discontinuous) nature, then at temperatures T c ~ 100 GeV there might 
have been a chance for developing the observed baryon asymmetry of the 
universe. The very active present day investigation of this question is based 
predominantly on a perturbative evaluation of the effective potential of the 
Higgs-field. 

The interest in non-perturbative studies stems from two sources - the 
generation of a magnetic mass as well as the breakdown of the perturba- 
tive treatment in the case of a weakly first order phase transition. A non- 
perturbatively generated screening mass of the magnetic part of the gauge 
field fluctuations is expected to weaken the discontinuous nature of the transi- 
tion. Eventually this could lead to a second order transition, if the magnetic 
screening mass becomes too large Thus the perturbative calculations 
[@, H lead to constraints on the maximal magnetic screening mass allowed 
for a first order transition and also lead to other predictions, like for instance 
the latent heat and surface tension, important for discussing the kinetics of 
the phase transformation. 

The validity of any perturbative treatment of the transition has been 
questioned for first order transitions where near T c the ratio T/m($ min ) be- 
comes large (weakly first order transition). In this case the fluctuations of the 
light modes (especially of the Goldstone modes) might require special con- 
siderations. Very recently various exploratory studies have been performed, 
in which it has been attempted to apply renormalization group theory of 
critical phenomena to this situation [|], |5], . However, one has to emphasize 
that techniques like the e-expansion were successful in finding universal char- 
acterisations, mainly critical exponents. Their success for calculating critical 
temperatures, order parameter discontinuities or surface tensions is by no 
means guaranteed. 

Under these circumstances numerical simulations seem to be an especially 
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valuable source of information. Exploratory investigations of the full, 4-d 
finite temperature theory met difficulties because of the large fluctuations 
due to the weakly coupled nature of the phenomena In this situation 
dimensional reduction might prove to be of important practical help @. P . 

In this work we study the influence of the finite temperature fluctuations 
of the Higgs-field on these parameters. In particular we will examine an 
effective 3-d theory for the Higgs-field, which has been obtained from the 
(3+l)-dimensional SU(2)-Higgs model in two steps: 

i) Dimensional reduction: One integrates over all non-static Matsubara 
fields at one-loop level. Since these modes are massive, no infrared 
sensitivity is expected, their perturbative integration seems to be well- 
founded. 

ii) Elimination of gauge degrees of freedom: The theory resulting from 
the first step is a 3-dimensional Gauge-Higgs model, where in addition, 
also an isovector field, the fourth component of the gauge fields, is 
present. In order to reduce the theory further one integrates over the 
magnetic gauge degrees of freedom and the isovector scalar. According 
to the improved perturbative treatments, the infrared stability can be 
ensured, if one includes into the result of the "naive" 3-dimensional 1- 
loop calculation a magnetic and an electric screening mass. Especially 
the first of them lacks, however, firm theoretical basis. 

We remark that the resulting effective scalar theory has been used re- 
cently also for the investigation of bubble nucleation in the electroweak phase 



transition [10]. 



In section 2 we first present the result of step i) for an SU(2) Higgs-model 
with N doublets: 

S = £dr J 

(1) 

(q = 1, N; fi = 1, ., 4; a = 1, 2, 3; Dfi = (d„ + igA*r a /2)<f>). 

In the limiting case N=oo step ii) can be performed exactly and a pure 
scalar action of very similar form to that actually studied in later sections is 
found. Also, some evidence will be presented from the perturbative analysis 
of the 3-d effective high temperature action of the N=l Higgs model that in 
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the coupling region corresponding to small Higgs mass values the pure scalar 
model might offer considerable insight into the physics of the electroweak 
phase transition. 

The discretisation of the effective model is described in section 3. The 
continuum limit of its mean field solution will be discussed in order to clarify 
the strategy of the non-perturbative investigations. Special attention will 
be paid to the uncertainties in taking the continuum limit, arising from the 
application of different renormalisation conditions. 

Section 4 presents a detailed discussion of the Monte Carlo simulation 
of the effective scalar model. Several physical quantities relevant to the 
phase transition will be stochastically evaluated and compared carefully to 
the information available from the literature. In the simulation the mass of 
the Higgs field was chosen approximately m# (T = 0) ~ 35 GeV. Conclusions 
will be drawn in section 5. 



2 Derivation of the Effective Scalar Action 

In this section we outline the derivation of the effective 3-d action for the 
SU(2) Higgs model with N scalar doublets on the 1-loop level. It represents a 
natural framework for the investigation of the finite temperature electroweak 



phase transition in an appropriate large N limit || [Hp. It turns out that 
the resulting 3-d action can be shown to be strictly equivalent to an effective 
pure scalar model. In this sense the model provides motivating background 
for the investigation of the scalar effective model in the case of a single scalar 
doublet (N=l). 



2.1 Effective Action in the Large N Limit 

Choosing the nonvanishing Higgs-field vacuum expectation value, <3>o> m the 
first scalar doublet, the computation naturally breaks up into a part identical 
with the calculation of the reduced action for the model with one doublet 
and the other part which consists of calculating the contribution from the 
appropriately separated group of N-l doublets. This statement is evident 
if one considers the quadratic part of the scalar action alone. The Higgs 
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doublets are parametrised as: 



Their potential is 0(4N)-symmetric and depends only on the length of the 
4N-component Higgs-field vector: 

Um sgs = \m^ q + ^A($+$ g ) 2 . (3) 

The quadratic action derived with eq.(Q) and eq.(^) has the form: 
i \ N 4 1 \ 1 \ 

q=2a=l A a^3 



Clearly, at 1-loop level the fluctuations of the fields with different q are 
independent, and this is true also for the gauge-scalar coupling part of the 
action. In the quadratic part of the latter the gauge fields couple only to the 
q=l doublet, where the static part of the Higgs-configuration was choosen to 
point to. The fluctuations of the other N-l doublets coincide with those of 
2(N-1) independent complex scalar fields on a static abelian A background, 
which is characterised in the Fourier-space uniformly by the following 2x2 
matrix: 



K 2 + m 2 + g -A 2 + igw n Ao 



(5) 



(K 2 = k 2 + u 2 ). The contribution from the coupled gauge-scalar fluctuations 



has been evaluated before f|, [L2| . We thus discuss here only the contribution 



from the additional scalar fields: 



AU = (N - 1) E / 7^3 ln ^ 2 + J^o + m 2 + ^l) 2 - g 2 u 2 n Al}. (6) 



n^O' 



(27T) 



The usual expansion in the argument of the logarithm leads to modification 
of the classical action, in which we are going to keep terms up to dim 4 combi- 
nations. After performing the frequency sum and the momentum integration 



4 



the following cut-off regularised corrections to the quartic $o — Ao-potential 
arise: 

AT , T 2 , m 2 \ , Xfh2/ A 2 AT 



AU = (N - l){f^(-AT + ^ + + _ 



2vr 2 



+ S - A,£) + ^ + ^Al$l + - D ^)} (7) 

with D = In ^+ const. Adding eq.(|7|) to the 1-do ublet4- gauge contribution, 
calculated in the thermal static gauge previously [fL2|| , the regularised full 
potential has the following expression: 

U Btatic {A ,$o] = ±m 2$ o+ 2l A$ o 

I 1^2^9.2 I 2N+1 x w A 2 AT , T 2 n r» 2 \ , m 2 / 3g 2 , (2AT+l)A x 1 

, 1 ^4/ (^+2) A 2 27g 4 , 3V _ n ( (N+2)\ 2 ?n£_ , 9A. 9 2 u 
24^0 1 12vr 2 64tt 2 ^ 16vr 2 12tt 2 64vr 2 ^ 16tt 2 >> 

+ l9 2 Aim + »^ + ^(5 _ 3A,)). (8) 
For the renormalised potential the "classical" extremum conditions are pre- 



scribed at the classical position of the minimum [[13 



dU(T -indep) cP[/(T - indep) 2 

^ = °' M = m n(T = 0), % = v . (9) 

For the separation of the T-independent part one has to rewrite the logarith- 
mic part (proportional to Do) of the regularised potential conveniently: 

_Ir^ m 2$2/9£ , 2N+1\\ , 1 ^4/ (^+2) A 2 2Jj£ , 9V^ , A 2 . 
2 L 167T 2 m 4 + 3 A ^ + 967r 2 ^0l 3 + 16 + 4 /J T 2 

= -64^E Q ^m 4 Q ($ )lnA!, (io) 

Q=T(ransversal),L(ongitudinal),H(iggs),G(oldstone). The most natural choice 
for mg($ ) is the following: 



n T 


= 6, 




n H 


= 1, 




n G 


= AN- 


1- 


n L 


= 3, 


ml 



m 2 + ^$ 2 , 



2 2 i ^f.2 

m G = m + 7T $ 0> 



A/ 
6~ 

27 \V2 



16 16(4tV-1) 



g 2 n (11) 
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with A' = A + 27 g 2 / (8N — 2). For large values of iV m 2 L is real, and the limit 
N — > oo gives 

n L = 3, m 2 = I^ 2 = m 2 ,, (N = oo) (12) 



On the right hand side of eq.fllOD one can then separate the cut-off dependence 
from the T-dependence by introducing appropriate normalisation scales, pos- 
sibly different for each value of Q. For the gauge coupling g 2 we choose the 
renormalization scale fi = T and absorb the logarithmic and finite corrections 
fully into its counterterm. This is consistent with the omission of the effect 
of the wave-function renormalisation in the present calculation (c.f. [|12]j ). 

After absorbing the T-independent infinities and some finite parts into the 
counterterms of m 2 , A and g 2 one finds the following explicit T-dependent 
potential fulfilling the conditions eq. @: 

U[A , J>„] = + i,A** + + ^±ai)(S - B) 

-T^- 2 *?!^ m |g* + A(4iV - 1) k ^ + fS(. + ^ + 2A] 



, 1 ^4r 27. 9 4 i ATgX _ 2(4iV-l) x/2 i ggf^d 
+ 15367r 2 ^0[ 2 111 4T 2 3 A 111 16T 2 

_ 6A 2 In ^ + A 2 ^ + f i7 4 + 1% 2 A] 

, 1 42/7V+4 2 T 2 , iVg 2 -m 2 _ AT+4„2 AT .\ i 3 4 , 4 4 , 1 2 /(2 < j ) 2 /io\ 

(($ 2 )t=o — Nvq). The large N form of this potential depends crucially on 
the way the different couplings scale with N: 

2 9 2 x A 2 m 2 
q = , A = — s-, m = . 14 

N a np N~y 

By inspecting the second line of eq.(Q) one realizes, following the argument 
of [|13|], that the stability of the symmetry breaking vacuum at T=0 requires 

^-<O(N ) or P-2a<0. (15) 
A 

This restriction excludes the proposition for the large N scaling law suggested 
in (/3 = 1, a = 1/3). We continue our analysis along the lines advocated 



in which corresponds to 

a = l, (3 = 2, 7 = 1. (16) 
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The length of the $o-field is assumed to be 0{N 1 ^ 2 ). We also note that 
the choice a = 1 for the scaling of the gauge coupling agrees with similar 
prescriptions used for SU(N C ) gauge theory coupled to Nf fermions in the 
Nf — > oo limit []I4 . 



Substituting these scaling relations into eq . (p~3|) , we obtain the following 
form of the leading large N potential: 

Ufa, $ ] = ^l\m\l - gSL) + (f + f )(g - g)] 
+ ^o(A + ^(l-lng|)) 

2^o(f^ 2 - fiAT) + ^ 2 ^«. (17) 



This result clearly shows that under the scaling conditions of eq.(|T6l) only 
the result of the gauge integration survives the large N limit. This conclusion 
agrees with the one obtained in Ref. |§ concerning the class of the leading 
Feynman graphs drawn in the full 4-d finite-T theory. Each term of the 
potential is O(N ), that is the same order of magnitude as that of the gauge 
part of the 3-d effective action. The latter one becomes the copy of three 
effectively Abelian gauge fields because g 2 now is scaling with an inverse 
power of N [IT]. The self-interaction of the 3-d isovector scalars A Q belongs 



also to the subleading terms of the potential, therefore its action becomes 
formally quadratic, too. With the (almost) usual 3-d-rescaling of the fields 
the leading N effective action is given by the following expression (the flavor 
index q has been omitted for convenience): 

| J <Px{U?[(-% + ~gl<P + <P)5 l3 + (1 - zWiP* + 2i9sVNU?J? 
U$[-d? + + § -^}U$ + <P + {-Nd? + m 2 ]0 + §§ (0+0) 2 } 

+ (3 — d) — counterterms (18) 

^2 _ ~2/i 27g| x ( 9g? s 2\ S \T \ _ \ I 27 3 4 (\ In tA\ 
771 - m l i MX^> + \— + —J 12' A 3 - A 3 + 128^11 m^j, 

g! = ~g 2 T, A 3 = AT, U? = A?/VT, U$ = A a /VT, = (19) 

The integration over Ui and Uq can be performed trivially. For instance in 
the Landau-gauge the non-regularised contributions to the effective 0-action 



with 
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reads: 

3Tr logi-df+gl^ + l Tr log(-5f + ^ + + f % J? D$J>, (20) 

where Df- denotes the propagator of the transversal gauge quanta on the 
scalar background. The last term is actually a correction to the kinetic 
term of the Higgs-fields. Taking into account that the fluctuations of <fi are 
0(l/y/~N), the magnitude of this term is also unity. After performing the 
3-d functional traces the 3-d mass-counterterms proportional to the gauge 
coupling exactly cancel and the following effective scalar theory is derived: 

S 3 M = J d 3 x{^+[-Ndf + m 2 ]0 + \NglJ?Dp* + |f (0+0) 2 

-U™a + ^) 3/2 - 1(1 + 4^) 3/2 - ( 21 ) 

In this equation we indicate, where could have appeared the magnetic screen- 
ing mass rh\, which to leading order in N is strictly 0. 

The structure of the effective theory eq. (|2T|) agrees with the basic struc- 
ture of the scalar theory to be investigated in later parts of this paper, up to 
the correction of the kinetic part which we shall omit from our quantitative 
investigation. 



2.2 The Effective Model for N=l 

For N=l one has to go back to the form of the potential term of the 3-d 
effective action given in eq.([U|). This form has been established in Ref. [fL^I , 
where also the 1-loop solution of the 3-d effective model has been extensively 
discussed. The authors suggest the existence of an interesting screening mass 
hierarchy at the temperature of the phase transition, whose characterisation 
depends crucially on the value of m#(T = 0). 

Namely, the effective screening mass squares of the A and A4 quanta 
given by the formulae 



^ = ^ 2 T 2 + ^ + ^(T = T c ), m^ = 9 ^l{T = T c ) (22) 



turn out to be a factor 4-6 larger than that of the scalar quanta for mn{T = 
0)=35 GeV. This ratio monotonically increases as one goes to smaller Higgs 
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mass values, its value is ~ 25 for m#(T = 0)=20 GeV. Therefore, at least in 
this region, the integration over these very massive degrees of freedom seems 
to be well-founded. 

After gaussian integration and after 3-d mass-renormalisation of A the 
following effective action for the Higgs-field is obtained: 

SsM = Jd 3 x[§^<W + |M|0+0 + f (0+0) 2 
-sf( 1 f i + «) 3/2 -^(^ + + 0) 3/2 ], (23) 

with the meaning of the notations given below (not to be mixed with previous 
meaning of some of them!): 



Ml = rn 2 + + -)T 2 - (24) 
* v 16 12 ; 2tt 2 ' 



m 2 = m 2 {l - 3^[(|^ 2 + A) In + A In M] _ ^(45y 2 + 20A + 

A = A-^ + ^ + f) 
-M9\l ^ gl - | In 0) + f In ^ + 3(¥ + f) 2 ^ (25) 

Eqs. (|23|) -(p5f) represent the version of the effective scalar model, which 
we are going to employ in the non-perturbative investigations of the next 
sections. Its general structure fully coincides with the result of the large N 
calculation, providing some extra argument for its relevance. 

For its discrete treatment it is most convenient to keep g 2 fixed and to 
think of the effective action as being defined by the couplings and A. 
After the performing the phase transition analysis one uses equations eq.(|2~4|) 
and eq. (p5|) for finding a posteriori the true couplings (or what is equivalent 
the Higgs mass). The deviation of the quantities with "hat" from "hatless" 
partners depends crucially on the logarithmic terms appearing in eq. (|24T ) 
and eq. (p5[) . With decreasing Higgs- mass they become increasingly impor- 
tant, making our conclusions for this region sensitive to the renormalisation 
scheme, implemented in the 1-loop calculations. For comparison we can con- 
sider a choice of the renormalisation scale in eq.(0), where the logarithmic 
terms are missing (/i = T). (The assumption of unchanged tree-level rela- 
tionship between the couplings and the W and Higgs masses is maintained 
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somewhat arbitrarily.) The results for T c and jf- will be seen in the next 
section to agree with their values found in the Linde scheme within 10-15% 
for rriniT = 0)=35 GeV when the mean field solution method is applied. 
The deviation sharply increases for smaller values of the Higgs mass. 

3 Lattice Formulation and Mean Field Anal- 
ysis of the Effective Scalar Model 

The discretisation of eq.(|23"D proceeds by introducing the dimensionless field 

i) = y/a(f> (26) 

with a being the lattice constant. Introducing also dimensionless masses and 
temperature via the relations: 

2m A a = y^fligQ, 2m D a = y^mgQ, 6 = aT (27) 

one finds 

_2ie^(2( TM e + ^xf' 2 + (lE<d + ^,) 3/2 )] (28) 



32vr 



with 



Here S(L 3 ) denotes the lattice regularized version of the linearly divergent 
integral J d 3 p/(8ir 3 p 2 ), giving rise to the linear divergence in eq. fl24|) . The 
perturbative values of the constants C and ^e are given by C = X and 
je — 10/3. We emphasize once more that the input parameters into the 
lattice calculation are k, X and g. 

For each fixed value of 9 the discrete theory, eq. fl2"8|) , will pass through a 
phase transition at a certain value k c (Q) or 

The continuum limit is defined by the limit — > of Z C (Q) providing a 
non-trivial value Zq. 



10 



When the value of Zq has been determined, eq.(29) is rewritten as a 
relation between the unknown physical value of T c and A: 

27f = Z °-2 ( 16^ + 12 } - (31) 



In eq.(31) rh is replaced by the expression given in eq.(25). The other relation 
is provided by the definition of A in eq. Q25D . From these two equations the 
physical values of the Higgs mass and of the critical temperature can be 
determined. 

In concrete terms, after introducing the dimensionless temperature (in 
units of the T=0 vacuum expectation value v$) the following two equations 
are to be solved (for instance, iteratively): 

A = A + / 2 (A,r c ), (32) 

with the functions 



/i(a, t c ) = t^(45^ 2 + 20A + *£) + 3^[(|^ 2 + A) In j£ + A In 



167T 2 V 16 1 4 1 3 

2 , \ . „ / 3„2 \ . „ _o„2 . 



+ s^ 4 (§ In £ - | In ^ ) + i In £ + 3(^ + |) 2 In %]■ (33) 

(Here the relation m 2 = — At> 2 /6 has been used again.) 

Having found A and r c the physical Higgs-mass and the critical tempera- 
ture are determined easily from the relations: 

rl = T l X = ^f. (34) 

The simplest illustration of a non-trivial continuum limit (Z ^ 0) is given 
by the mean field analysis of the system eq. fl28|) . The mean- field expression 
of the free energy density as a function of the mean field order parameter s 
is 
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The condition for the existence of a non-trivial (sq 7^ 0) minimum is written 

as 



Introducing the scaled free energy F[s]/G 3 and the scaled order parameter 
s 2 = s 2 /Q, one can eliminate the common term (j- — rf)^ from the degen- 
eracy condition for the free energy minima and determine the position of the 
degenerate symmetry breaking minimum Sq from the equation 

-i4 + ¥#[(7m + 4) l/2 + Kie + 4) 1/2 \ 

-£[2( lM + 4? /2 + (IE + 4) 3/2 } + £(2j¥ + if) = 0. (37) 
The quantity Z c>m f is calculated as 



We note that the right hand side is independent of 9 and therefore Z , m / = 
Z C}Tn f. In Fig. 1) we show the mean field estimate of the quantity Z for 



7s = 0, 1, 2 (see eq.(f4T|) for the meaning of jb) as a function of rhu- The 
rapid rise for small % leads uniformly (independent of 7#) to a "critical" 
A at which the right hand side of eq. (|3lD vanishes. In the approximation 
m? = m 2 and A = A this corresponds to a divergent critical temperature. 
However, because of the uncertainties due to the application of different 
renormalization conditions this phenomenon might be an artefact of the 
above approximation. The data point appearing in Fig. 1) is the result of 
our numerical simulation. The analysis leading to it is described in section 
4. 

The order parameter discontinuity characterising the phase transition can 
be expressed with help of so as 

$(T C ) = s T c . (39) 

A very contentful characterisation of the first order transition is given by 
the surface tension, a, between coexisting ordered and disordered regions. 
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Exploiting the scaling behavior of the mean field solution one finds for it in 
the thin wall approximation: 



a mf = ^ = f°°ds[2f(s 

c 



f(s) = = Z o, mf s 2 + ^ - £[2( lM + s 2 f' 2 + { lE + 

+ £(2l% 2 + lT)- (40) 

In Fig. 2) we show the true Higgs mass square as a function of rh H = 
A^q/3. The trivial linear relation goes over into a strongly non-linear func- 
tional form only for Higgs mass values less than 25 GeV. For larger values 
niu an d run are seen to agree within 10%. In Fig. 3) T c is shown as a func- 
tion of mil and m#. Below mu =25 GeV T c is rising with decreasing Higgs 
mass, reflecting the approach to zero of the denominator in the expression of 
r c in eq.(|32]). 

The solid curve in this figure shows the result of the approximation 
mu — Ah, A = A. This approximation corresponds to the choice of nor- 
malisation scale fi = T in the interpretation of eq . (|T0|) . It leads to the 
definition of the continuum limit directly from eq.(|3~I|) avoiding the compli- 
cated procedure of solving eq.(p2[) and eq.([33|). This approximation for a 
Higgs mass near m#=15 GeV would indicate infinite T c , which would hint to 
the non-restorability of the broken gauge symmetry for Higgs masses smaller 
than this limiting value. In principle, the non-restoration of the symmetry 
even for infinite temperature is not excluded and there are known examples 
of such behaviour due to quantum corrections to the effective potential [|15l 



These are the logarithmic terms of the expression of A and m 2 , which counter- 
balance the increasing tendency of T c , when the Linde-type renormalisation 
condition is being used. However, large logarithmic corrections due to the 
essential deviation of the relevant range of the temperature from the chosen 
renormalisation scale make the 1-loop approximation unreliable. The differ- 
ence in the two curves of Fig. 3) for mu < 30 GeV indicates the importance 
of higher loop contributions in the equations governing the continuum limit. 

In Fig. 4) the dependence of the order parameter discontinuities on the 
strength of magnetic screening is shown as a function of m#. For non- 
vanishing magnetic mass the first order transition ends in a tricritical point 
(TCP). Following we parametrize the amount of magnetic screening in 
proportion to the value obtained from a selfconsistent Dyson- Schwinger equa- 
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tion: 

g 2 T 

m A = 1b — - (41) 

The mean field analysis yields then for the end point values of the Higgs- 

masses: 

7fl = 1 , m H (TCP) = 84.5 GeV, 

7b = 2 , m H (TCP) = 61.7 GeV (42) 

in good agreement with the continuum perturbative estimates. In the inser- 
tion of Fig. 4) the uncertainty in the order parameter discontinuity resulting 
from the use different renormalisation schemes is illustrated. The application 
of the Linde-conditions leads for small Higgs masses (dotted line) to harder 
transitions. We can conclude, that in the region m^j < 30-35 GeV no reliable 
continuum statement can be made on the basis of the 1-loop calculation. 

Finally we display in Fig. 5) as function of % the mean field result for 
the interface tension eq. (|40|) . 



4 Numerical Simulation and Results 



We have simulated the lattice model eq . (p8|) on hypercubic lattices with 



linear extent L and periodic boundary conditions. The electric and magnetic 



screening mass parameters were choosen to be 7^ = 10/3 and 7^ = 1, eq. (|41"|) . 
At values for the W-mass my/ = 80 GeV and for the zero temperature 
vacuum expectation value of the Higgs field vq = 246 GeV the gauge coupling 
constant g of eq.([2"8"D is determined by the tree level relation eq.(^) applied 
at T = 0. These choices first leave us with 3 free parameters of the theory 
A, B and the hopping parameter k. Choosing furthermore a value for the bare 



quartic coupling A or correspondingly via eqs.( p2] , p3| , p4D a physical value for 
the Higgs mass m#, we are left in our simulation with a two parameter 
theory. In the k-0 plane of couplings we expect the symmetry restoration 
phase transition whose properties we would like to study. In detail we shall 
discuss how to extract T c , the order parameter discontinuity, $(T C ), the latent 
heat £ and the interface tension a from the raw data. 
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4.1 Lattice Observables and the Continuum Limit 



Both bare parameters k and 6 serve the purpose of constructing the contin- 
uum limit. In particular we will be interested in the study of various quanti- 
ties measured at T c . We thus will tune the hopping parameter to the critical 
line k c (G) and consider the limit —>■ 0, which at the same time removes 
the lattice cutoff a. In this limit the theory turns into the 3-dimensional 
Gaussian model. From the 0-dependence of the critical hopping parameter, 
we can determine the continuum limit of the quantity Z C (Q), which has been 
defined in eq.(^). With the help of eqs.(^,^,0) this fixes the critical tem- 
perature in units of the Higgs mass. We furthermore study various operators, 
which allow the extraction of physically relevant quantities such as the order 
parameter discontinuity at T c , the surface tension and the latent heat. In 
particular we consider the following observables 



whose distributions P(Oi), expectation values (Oi) and diagonal fluctuations 
C(Oi) = L 3 (0 2 — (Oi) 2 ) for i = 1,2,3 are determined in the simulation. 0\ 
corresponds to the average length square of the Higgs field, O2 to a part 
of the kinetic term of the action, while O3 corresponds to the average field. 
They can be related to the corresponding continuum observables through 
relations discussed in the previous sections. Here we note that the lattice 
fields ip and the continuum fields $ differ by a factor of VT. In the limit 
— > we thus expect appropriately scaled dimensionless ratios to reach 
constant values. For instance we examine the ratio Qiyi, which in the limit 
0^0 converges to the continuum result $(T C )/T C . Of course we have 
to investigate here carefully possible scaling violations due to finite lattice 
spacing a, as well as finite volume effects due to the finite extent L 3 of the 
lattices. These problems will be discussed in detail in connection with the 
numerical analysis of the various observables. 

In general we expect that with decreasing values of (decreasing lattice 
spacing) it will be increasingly difficult investigate the properties of finite 
temperature phase transition as any signal for a possible first order phase 
transition will become small in units of the lattice spacing. Moreover, we 
have to understand the possible range of values of the Higgs mass, for which 
we can study the phase structure of our model. For large values of the Higgs 




x+e ) 




(43) 
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mass the transition is expected to be at best weakly first order, which then 
requires large lattices to resolve any discontinuity in physical observables. In 
order to explore the parameter range, in which we can perform statistically 
significant studies of our model we have studied thermal cycles across the 
phase transition line k c (Q) varying also the Higgs mass parameter. 

In Fig. 6) and Fig. 7) we display the results of thermal cycles for the 
order parameter (0 3 ) on 18 3 lattices with a statistics of 4000 sweeps per data 
point. Fig. 6) displays results at fixed values of = 5. The parameter rhn 
ranges from m# = 35 GeV to rhn = 47 GeV in this figure. The broad and 
pronounced hysteresis at small values of rhu rapidly shrinks with increasing 
wiff and disappears at rhn = 47 GeV. We have included into both figures the 
mean field result for the critical point (dashed vertical lines) and the mean 
field value for the order parameter discontinuity (solid circles and height of 
the dashed lines). While the locations of the phase transitions roughly agree 
with the centers of the hystereses, Fig. 6) indicates a more rapid decrease of 
the order parameter jump as compared to the mean field result. At rhn = 47 
GeV and = 5 it is presumably notoriously difficult to determine infinite 
volume values of possible gaps in operators, or to determine whether the 
phase tranition is discontinuous at all. This is the main reason why we have 
decided to choose the value rhn = 35 GeV for a detailed and high precision 
study of the symmetry restoration phase transition. This value of rhn = 35 
GeV is still large enough for not being plagued too much in the continuum 
interpretation by the ambiguities of different renormalization schemes. It 
corresponds to a physical value of the Higgs mass of m H = 37.16 GeV, see 
Fig. 2). 

Fig. 7) displays at the selected value of m# = 35 GeV the variation of the 
hysteresis of the order parameter O3 as a function of 0. As expected the first 
order signal weakens with decreasing 0. Comparing the mean field values 
of the order parameter jump with the MC data we again note a more rapid 
weakening of the first order signal. From these data and from additional 
MC-runs we expect that the interesting region in which we might be able to 
resolve possible first order signals on the given sized lattices extends down 
to 0- values as low as = 1. The typical lattices sizes, which we use in 
our simulation range from 8 3 up to 18 3 . Typical statistics for simulations 
on any of the lattice sizes is about 10 6 sweeps. To overcome the problem of 
large tunneling times at the first order phase transition, we have used the 
multicanonical ensemble approach We see later that even at = 3 
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the phase transition appears to be rather week, e.g. double peak structures 
in the probability distributions P(Oi) at the transition point are not very 
pronounced. 

4.2 The critical Temperature T c 

The perturbative treatment of the effective model predicts a certain shape 
of the phase transition line k c (Q), provided Z C (Q) of eq. (|30|) is known. How- 
ever the mean field treatment of the theory predicts that Z C (Q) = Z 0jTn f 
is just a constant independent of O. We therefore might expect that in a 
first approximation also the fully fluctuating theory can be described with a 
constant Z C (Q) = Z ^ mc for small values of 6. At rhu = 35 GeV we estimate 
the critical hopping parameter from thermal cycles. These are displayed in 
Fig. 8) for less than 1. We observe a clear bending of the phase tran- 
sition line towards smaller K c -values at O-values of about unity. Thus the 
numerical data predict a nonzero and positive value for Z 0jmc . A fit to the 
shape of the phase transition line according to eq.(|30|) then results into a 
Monte Carlo determination of the quantity Z rnc = 0.0187(16), which is to 
be compared with the mean field value of Zo,m/ = 0.0113 at same values of 
couplings. Noticeably our data are consistent with the presence of a linear 
©-dependent term in eq. (|30|) with C = 0.93(16) x A , as it is predicted by 
perturbation theory (the continuous curve in the figure). The mean field 
phase transition line (dashed curve) deviates from the MC data, which at 
small G slightly overshoot the k- value k = |. At the given value of Z ^ mc we 
may now determine the critical temperature for both of our renormalization 
schemes, namely in the \i = T scheme (without logarithmic corrections), as 
well as in the Linde-scheme, which incorporates them. The result is 

T c = 114.9(36) GeV (jj, = T) 

T c = 114.3(30) GeV Linde scheme. (44) 

Both numbers are displayed in Fig. 3), the triangle corresponds to the \i = T- 
scheme and is plotted at its value of rh H while the circle corresponds to 
the renormalization scheme with logarithmic corrections plotted at the the 
physical Higgs mass value m#. It is noticeable that the difference between 
the schemes is very small. These Monte Carlo determinations of T c should 
be compared with the mean field result T c = 99.6 GeV. 



17 



4.3 Critical Hopping Parameters 

For the value of rhu = 35 GeV, we have performed a detailed finite size 
scaling analysis of the symmetry restoration phase transition at three values 
of G, namely = 3, 4 and 5. In Fig. 9) we compare the results of a thermal 
cycle on an 18 3 lattice for the operator (0\), with results of multicanonical 
simulations for the same operator on 8 3 , 10 3 and 12 3 lattices. The mean 
field position of the critical point, as well as its true infinite volume value, 
deduced from finite size scaling analysis, are indicated by the vertical dashed 
lines in the figure. Defining pseudocritical points K max (L) from the peak 
positions of the maxima of the diagonal fluctuation C(Oi), see Fig. 10), the 
pseudocritical ^-values are extrapolated to infinite volume via 

Kmax(L) = K c + — . (45) 

Fig. 11) shows at = 4 the finite size scaling analysis leading to the infi- 
nite volume critical point. It is clear from these figures that multicanonical 
ensemble simulations combined with finite size scaling theory allow a precise 
determination of the location of the phase transition or discontinuities. They 
are in fact superior to standard methods e.g., thermal cycles, which even on 
a large 18 3 lattice cannot produce numbers with comparable precision, see 
Fig. 9). Table 1) then contains for the selected 0- values a comparison of the 
mean field critical hopping parameters to the Monte Carlo results. They are 
very close and their difference is O(10~ A ). As one might have expected the 
location of the phase transition at large values of is well described by mean 
field behavior. This has to be contrasted with the observed behavior of k c (Q) 
at values of < 1, where the mean field approximation is not adequate. 

4.4 Determination of Discontinuities 

We obtain the infinite volume discontinuities in thermodynamic quantities by 
analyzing the volume dependence of probability distribution functions -P(Oj) 
of operators Oi defined in eq . (f43[) . In the vicinity of the transition they 
develop double peak structures. We determine the values of the couplings, 
where in the finite volume system the two peaks in P{Oi) are of equal height. 
This defines e.g., at fixed value of 0, a value of the hopping parameter n e h(L). 
Fits to the maxima of the distribution functions then yield finite volume 
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Table 1: The critical hopping parameters K C (Q), to the left: mean field 
results, to the right: MC results. 






Kc(Q)m/ 


^c(0)mc 


3 


0.161188 


0.161112(4) 


4 


0.157170 


0.157065(5) 


5 


0.152283 


0.152170(2) 



estimators for the values of operators in the metastable states, denoted by 
{Oi)sB{L) (symmetry broken) and (Oi)s(L) (symmetric) for i = 1,2,3. The 
discontinuities A(0$)(L) = (Oi)sB{L) — (Oi)s(L) are then extrapolated to 
infinite volume using standard finite size scaling arguments ]T7|. Their large 
volume L-dependence is expected to scale with the volume, 

A(Oi){L) = A(O t ) + i = 1,2,3. (46) 

In Fig. 12) we display on 10 3 lattices probability distributions of the operator 
0\ i.e., the average length square of the scalar field at n e h(L) for values 
= 3,4 and 5. As can be seen one finds clear double peak distribution 
functions, which can actually be observed for all of our operators Oj % — 
1, 2, 3. Remarkably, however, the value of the minimum of the distribution 
functions rapidly rises with decreasing 0. At = 3 the suppression of mixed 
phase states is hardly noticeable, which makes the phase transition a rather 
weak first order phase transition already. 

Fig. 13) shows the finite size scaling analysis of (0 3 )sb(L) corresponding 
to the order parameter jump in the symmetry broken phase. A fit to the data 
with the form eq.(|46|) readily gives the infinite volume discontinuity A(0 3 ). 
Fig. 14) shows the same kind of analysis for the operator Oi, corresponding 
to the average length square of the scalar field. In this case (Oi)s(L) also 
has to be considered for the calculation of the corresponding discontinuity. 
Carrying out the analysis for all our operators we arrive in Table 2) at the 
infinite volume values for the discontinuities for all considered operators and 
0-values. They can be compared with the mean field order parameter jump 
which appears also in Table 2). In Table 2) we present the discontinuities 
in 0\ and 2 scaled by the square of the discontinuity in 3 , the order 
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parameter. They are consistent with a value of 1. Thus actually the various 
discontinuities are not independent at the transition, but trivially connected. 
This is exactly what one expects from a mean field type of behavior at the 
transition. 

The figures also exhibit some unusual finite size scaling behavior e.g., it 
can be noted that with increasing lattice size the finite volume estimators 
of the discontinuities increase. Many of the first order phase transitions in 
statistical physics and lattice gauge theories exhibit an opposite behavior, 
that is with increasing lattice size the values of the obtained discontinuities 
decrease. Figs. 13) and 14) also indicate that finite volume corrections 
become rather large in size at the smaller values of G (see the different slopes 
in the same figures). This is especially true for the symmetry broken state 
and we attribute this property to the vicinity of the Gaussian fixed point. 



Table 2: Values of the discontinuities. 






Smf 


A(0 3 ) 


A(O a > 2 


A(0 2 > 
A(O a > 2 


1 


1.312 








3 


2.285 


1.603(50) 


0.955(70) 


0.946(69) 


4 


2.639 


2.036(05) 


0.998(07) 


0.999(07) 


5 


2.951 


2.458(05) 


0.991(06) 


1.000(07) 



4.5 The Order Parameter Jump 

In Fig 15) we show as a function of the order parameter jump A(Os) scaled 
by the mean field result. The actual magnitude of the order parameter jump 
comes out smaller, than predicted by the mean field analysis. In addition 
we also observe a mild deviation from scaling and it appears as if possible 
scaling deviations might be parametrized by a correction linear in 0. Fig. 15) 
contains a corresponding fit (dotted line in the figure), which predicts that the 
order parameter jump in units of the mean field order parameter discontinuity 
might in the continuum limit — > be as low as one half. However, the 
theoretical status of scaling corrections in 3-dimensions is unclear [|18|] and in 



our context we might regard such an extrapolation as a lower bound to the 



20 



true value of the order parameter jump. Further simulations closer to the 
continuum are therefore required. They may show a bending of the order 
parameter discontinuity to somewhat larger values than half of the mean field 
result, but will be significantly lower than the mean field prediction. Such 
a result is consistent with the variational upper bound nature of the mean 
field approximation. Employing our extrapolation we currently estimate the 
order parameter discontinuity to be 

^£p) = 0.68(4). (47) 

-L c 

Using the T c determination we obtain $(T C ) = 78(5) GeV at a physical Higgs 
mass value of m H = 37.16 GeV. 

4.6 The Latent Heat 

The volume normalized latent heat in units of T c 4 , is obtained from the 
expectation value of the internal energy ^ 

<^> = <^W)> (48) 

by calculating its difference in the symmetric and the symmetry broken state 
at T c . 

= (^j) Is -(^4) \sb (49) 

c c c 

Hereby f3 denotes the inverse temperature (3 — ^ and S denotes the action 
of the 3 + 1 dimensional theory. We note that our convention in defining 
the latent heat is opposite to what we have used for the order parameter 
discontinuities. This makes the latent heat positive. Incorporating the full 
temperature dependence of all of our physical couplings we obtain for the 
internal energy in terms of the 3-dimensional reduced theory the expression 

lU_\ _ I S 3D ,iat 1 v 1 A _ m H (In 2 -I- \\ tilts 

\ T 4/ — \ G 3 -r l^xl { 8 -t- 6 167r 2 T 2l 4 y -r a) @ 



_fM£ _i_ dl _i_ jA g y wfe ^i 

32tt 2 ^ 16 ^ 3 ^ 4 > v O / 
^)l + |(f + 

+S(1f) 2 ((^f ) 2 + + s(t + ^)']}> ( 5 °) 



+i4[(( 2 if ) 2 + + Kf + ^) § ] 
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We note here that the latent heat calculation requires the calculation of 
discontinuities of more complicated operators e.g., the calculation of discon- 
tinuities in J^xO'x^x) 2 and discontinuities in certain functions of the fields. 
Since we did not measure all required operators directly, we assume in the 
numerical analysis that corresponding discontinuities show a mean field like 
behavior e.g., the discontinuity in the operator E^l^i) 2 ) for instance, is 
replaced by the discontinuity in the operator 3 to the fourth power. Our 
data support such an assumption, at least, for the restricted set of operators 
at the large values of 6, considered in the present investigation, see Table 2) 
and the related discussion. 

In Fig. 16) we display the discontinuity of the volume normalized 3- 
dimensional action eq.(^) scaled by corresponding B-powers at the critical 
point (circles) as a function of 9. It is consistent with being zero and results 
from a cancellation of positive and negative terms, which can be attributed 
to the kinetic term and the potential term in the action. As the action gap 
itself is related to the pressure, which is continous even at the first order 
phase transition, this is an expected behavior and merely is a check on the 
consistency of our analysis. 

The latent heat per unit volume ^ also appears in Fig. 16). The triangles 
denote the results of our Monte Carlo simulation, while the crosses (constant 
with 0) correspond to the mean field analysis of the theory. Again we observe 
a significant weakening of the phase transition as compared to the mean field 
treatment. We also observe scaling deviations i.e., as a function of G the 
latent heat values have not yet settled to their continuum value. Here we 
quote the value of the latent heat obtained from simulations at = 3: 

£ = 0-122(8). (51) 

c 

This value again is less than half of the mean field result ^ = 0.262 and 
should be close to the continuum result provided further scaling deviations 
turn out not to be too large. 



4.7 The Interface Tension 

Following Binder [BJ the interface tension inbetween symmetric and sym- 
metry broken states of the theory can be obtained from the shape of the 
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probability distribution functions at n e h. Denoting the value of the maxima 
of the probability distribution functions P max (Oi), states inbetween the two 
maxima correspond to mixed phases, which on the largest lattices form an 
interface with a cross sectional area L 2 . Because of the periodic boundary 
conditions there are actually two interfaces in the system and the total in- 
terfacial area thus is 2L 2 . Denoting the value of the probability distribution 
function at its minimum P mm (Oi), the interface tension can be defined by 
the limit 

The quantity oj ot is connected to the continuum interface tension a via 
&iat — w® 2 - It * s however expected that the extrapolation to infinite volume 
is more subtle, than in the case of discontinuities in the order parameter. In- 
terfaces in finite boxes show sizable fluctuations controlled by their stiffness 
and finite volume corrections to <Ji at may have a complicated L-dependent 
analytical form, as is demonstrated by numerical simulations in spin models 
In Fig. 17) and Fig. 18) we exhibit our interface tension analysis. 
Fig. 17) displays the quantity ln( p mi „|^j ) for three O-values as a function 
of the expected interfacial cross section L 2 for all considered lattice sizes. 
In particular for = 5 we observe a very fast rise with L 2 , indicating that 
mixed phase configurations are suppressed by many order of magnitudes in 
the path integral and that the interface tension thus is large. The data are 
also roughly consistent with a linear increase of the considered quantity with 
the interfacial cross section, though some curvature to larger slopes is in- 
dicated. Since the status of possible finite volume correction terms to our 
measured quantities is unclear, we may regard the straight line fits of Fig. 17) 
as estimates for the expected surface tensions values. They are displayed in 
Fig. 18) as functions of in units of the corresponding mean field result ob- 
tained in the thin wall approximation, eq.(|40|). They show indeed a dramatic 
difference. In the mean field calculation we expect at the given couplings a 
behavior according to a\ a t = 0.0242 2 . Definitely large scaling deviations 
are observed, while at the same time e.g., at = 3, our estimates are about 
37 times smaller than the corresponding mean field calculation result. Thus 
it is indicated here, that the interface tension is subject to large corrections 
due to the fluctuations of the scalar field. Turning our interface calculation 
into physical numbers the data point at = 3 predicts an interface tension 
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of as low as about 960 GeV 3 , which in the current context may be viewed 
as an order of magnitude estimate. The numerical situation is not totally 
satisfying and in the future one should also explore other possibilities for 
interface tension determinations, like simulations on asymmetric lattices and 
studies of tunneling masses on such geometries. 



5 Conclusion 

In this paper a pure scalar model has been proposed for the description 
of the finite temperature electroweak phase transition (EWPT). This effec- 
tive model was obtained from the full 3+1 dimensional field theory through 
integration over the non-static Matsubara modes, realised in the 1-loop ap- 
proximation. The remaining 3-dimensional gauge degrees of freedom and 
the effective adjoint Higgs-field were eliminated in a subsequent Gaussian 
integration step. 

Intuitive justification for such scenario is provided by the thermal mass 
hierarchy observed in the perturbative treatment of EWPT for moderately 
small Higgs masses. Also, we have demonstrated that in an appropriately 
defined large N limit the second integration step is exact. 

The resulting model is a 3-dimensional 0(4)-invariant ferromagnet with 
cubic and quartic potential. 

Its discretisation and the proper procedure for taking the continuum limit 
of its phase transition have been thoroughly explained. Especially, the sen- 
sitivity of the T c dependence on m#(T = 0) stemming from the application 
of different renormalisation conditions in the continuum has been discussed 
in detail. 

The mean field solution of the lattice system and its continuum limit 
has reproduced in all aspects the results of the improved 1-loop perturbative 
treatments of EWPT [2,3] (influence of magnetic screening, order parameter 
discontinuity, surface tension, etc.). Also it called attention to the dramatic 
sensitivity of T c to the renormalisation conditions for m#(T = 0) < 30 GeV. 
This circumstance has prevented the discussion of the possible existence of a 
lower limiting Higgs mass value below which the gauge symmetry could not 
be restored at any finite temperature. 
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The mean field analysis served for orienting the numerical calculations, 
performed at m#(T = 0) ~ 35 GeV. This value has been selected : i) to have 
a strong signal of a first order transition, ii) to be free of the uncertainties 
introduced by the renormalisation prescriptions , iii) to be able to compare 
the results of a careful numerical study with other computer studies of EWPT 
based on different effective 3-d models. 

The numerical data consistently exhibit at the considered value of the 
Higgs mass a weaker first order phase transition as compared to the zeroth 
order mean field approximation of the theory. Discontinuities of the theory 
like the order parameter jump and the latent heat turn out to be somewhat 
smaller than corresponding mean field results. The critical temperature T c = 
114 GeV appears larger than mean field, again consistent with a weakening 
of the first order phase transition. The interface tension is strongly affected 
by the fluctuations of the scalar field and its value will be much lower than 
corresponding mean field predictions. By themselves these findings may not 
be regarded as surprises as naively one expects the fluctuations of the Higgs 
field to work in the observed direction. Similar experience has been gained 
from numerical estimations of the interface tension of the SU(3) pure gluon 
theory. 

Comparing our results with the results of another MC simulation making 
use of another variant of 3-dimensional effective models of the electroweak 
phase transition at comparable value of the Higgs mass, we make, however, 
an interesting observation. While the authors of |§ keep the gauge fields 
in the 3-dimensional Gauge-Higgs model as dynamical degrees of freedom 
in their simulation, their findings indicate a strengthening of the first order 
phase transition, when compared to a one-loop perturbative analysis of their 
model, and also, when compared to the results of our simulation. They 
do quote e.g., a lower value of T c = 85 GeV at a value of m# = 35 GeV 
or a value of the order parameter discontinuity of about 170 GeV, while 
our value <&(T C ) = 78(5) GeV is about half as large. Thus it appears as if 
the non-perturbative dynamics of gauge fields in the dimensionally reduced 
model works in a direction opposite to the effect of the fluctuations of the 
Higgs field, and eventually dominates the physical characteristics of the phase 
transition. 
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